Stripe as an effective one-dimensional band of composite excitations 
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The microscopic structure of a charge stripe in an antiferromagnetic insulator is studied within the 
t- Jz model using analytical and numerical approaches. We demonstrate that a stripe in an antiferro- 
magnet should be viewed as a system of composite holon-spin-polaron excitations condensed at the 
self-induced antiphase domain wall (ADW) of the antiferromagnetic spins. The properties of such 
excitations are studied in detail with numerical and analytical results for various quantities being 
in very close agreement. A picture of the stripe as an effective one- dimensional (ID) band of such 
excitations is also in very good agreement with numerical data. These results emphasize the primary 
role of kinetic energy in favoring the stripe as a ground state. A comparative analysis suggests the 
effect of pairing and collective meandering on the energetics of the stripe formation to be secondary. 
The implications of this microscopic picture of fermions bound to the ID antiferromagnetic ADW 
for the effective theories of the stripe phase in the cuprates are discussed. 
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I. INTRODUCTION 

Strongly-correlated models of the Cu02 planes of high- 
Tc superconductors continue to attract much attention 
due to the belief that most of the physics in the cuprates 
is governed by strongly interacting, purely electronic 
degrees of freedomtl. Microscopic studies of Hubbard 
and t-J models have been successful in expLainiug the 
d-wave character of the pairing mechanismodij^El and 
other experimental results, such as narrow low-energy 
bands in the anele-resolved photoemission for the un- 
doped systems.l3iZl'El More recent interest in these models 
has been boosted by the discovery of stripes, cuj-Sipin and 
charge inhomogeneities, in high-Tc materials.oEj Gener- 
ally, the strongly correlated t-J and Hubbard models in 
the non-perturbative regime {J t or U ^ t) are dif- 
ficult to approach analytically, although some advances 
have been achieved in solving them. Because of this, nu- 
merical methods have been responsible for murh. af the, 
progress in the understanding of these models.EJI13o'Efl 
Moreover, such numerical studies have become a very 
important test of the ability of theoretical approaches 
to describe the stripe and other low-energy phasealj in 
cuprates. 

While the striped phase was anticipfaied from mean- 
field solutions of the Hubbard modelll3, probably the 
most convincing evidence of stripe-like ground states has 
been provided by Density Matrix Renormalization Group 
(DMRG) studies of the t-J model in large clusters in the 
range of parameters relevant to real systems.EZI However, 
some other numerical approaches raise the question that 
the stripeSr-seen in DMRG might be the result of finite- 
size effectstS. Another aspect of the problem is that nu- 



merical studies alone do not directly answer questions on 
the origin of stripes. Ideally, one would wish for a the- 
ory which would closely agree with the numerical data 
on all essential aspects, thus providing a definite physi- 
cal answer on how the stripes are created and what are 
the excitations around this state. The search for such a 
case has motivated our preserttjwork and that is the way 
we unify our approaches here £3 

In this work we attempt to integrate some of the ear- 
lier ideas on the t-J model physics with the newer trends 
and phenomenology which have appeared due to stripes. 
We approach the problem using a comparative study of 
the stripe in an antiferromagnetic insulator by DMRG 
and an analytical technique, within the framework of the 
t-Jz model. Our numerical study utilizes DMRG in large 
Lx X Ly clusters of up to 11x8 sites, using various bound- 
ary conditions. The analytical method is a scl£-«orisistent 
Green's function technique developed earlierE3Eil, which 
accounts for the retraceable-path motion of the holes 
away from stripe. We demonstrate that the stripe in 
an antiferromagnet (AF) should be viewed as a system 
of composite holon-spin-polaron excitations condensed at 
the self-induced antiphase domain wall ADW. 

The t-J model has long been seen as a natural 
model for th£-|description of the charges and spins in 
a doped AFE3. The single- and two-hole problems 
within the model have been studied extensively using 
different analytii 
in small clustersi 

attempts to generalize the conclusions 
ies on the natutg. of the many-hole ground-state have 
also been madajSoEj. Very good agreement be- 
tween the numerical and analytical studies for these 



2aJ2'n,E8lG9 



roaches 
1^. Some 
of these stud- 
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problems. b.as. be^n achieved within the spin-polaron 
paradigmEJcjO'E^. This essentially quasiparticle picture 
describes the single-hole excitation as a hole strongly 
dressed by the "string" -like spin excitations. It was found 
that two such quasinar.tinkjs t.ond .tn.£Qi-m a bound-state 
of d-wave symmetry.yy'y'BEyO'E3 The generic rea- 
son for the absence of the s-wave pairing is the magnon- 
mediated exchaHftQ— which generates a repulsion in the 
s-wave channel .□u'EH—Lel the t-J model, vertex correc- 
tions are suppressecoo and such a repulsion is strong. 
Then the attraction in higher-order .hfirmonics leads to 
the bound states of higher symmetryc^l. Attempts to in- 
tegrate out the spin background and to reformulate the 
t-J model as an effective model for quasiparticles with 
a narrow band (~ 2 J) and an interaction repulsive in 
the s-wave chaniiel and attractive in the d-wave channel 
have been madeJj'OoEEl assuming the antiferromagnetic 
correlation length to be the largest scale in the problem. 
Phase separation in such a model in a physical range 
t 3> J seems to be unlikely since the pair-pair interaction 
should also be repulsive, which agrees with the numerical 
datac3. The ground state in this model would be a dilute 
"gas" of d-wave spin-polaron pairs. 

Such a generalization of the spin-polaron picture to 
a finite concentration of holes relies on the assump- 
tion that the antiferromagnetic background remains un- 
changed. However, it is well known that the "feedback" 
effect of holes on the antiferromagnetic background is 
important. Aside from Hartree-Fock treatments of the 
Hubbard modelE£l which showed stripe-like domain wall 
solutions, other studies of the t-J model in the low- 
doping regime hayei indicated instabilities of the antifer- 
romagnetic order.clThese instabilities^ were thought ta 
lead towards spiraCj, stripe-like spiralo, or spin-liquido 
states. Earlier numerical works in the small t-J clus- 
ters, Ref. [Q, have demonstrated stripes in the ground 
state which were also domain walls in the IXfeJ AF. With 
the mounting evidence-from experimentaa'ta and from 
DMRG numerical dataO the idea of topological dopingca 
has flourished. The spontaneously created ADW's have 
been widely considered as the topol og | iC | al . alt ernatives to 
the homogeneous Neel background.oBOS 

Thus, the many-hole ground state has turned out to 
be very different from that for a few holes. In order to 
understand the nature of the charge excitations in this 
phase one needs to reconsider the single-particle pmblem 
around this ground state with different topology.u The 
one-dimensional (ID) character of the charge stripes has 
led to a number of attempts to generalize the physics of 
strictly ID systems, where the (pajdlations are holons and 
spinous, to higher dimensions .OLJ On the other hand, 
there is a growing understanding that the stripes are the 
outcome of the same tendeocies which are seen already 
for the single-hole problemE3, and that the charge exci- 
tations in the stripe phase may still have lots in common 
with the spin polarons.EllEll 

Note that topological doping generally refers to the in- 
troduction of dopants into the topological defects of a 



field theory, which are the field configurations that inter- 
polate between different vacua of the problem. For prob- 
lems with high symmetry, such as the SU (2) Heisenberg 
model, topological defects are continuous distortions of 
the order parameterEa The ADW's are the topological 
defects for systems with lower Z2 symmetry, such as the 
Ising or the axiisotropic t-J^ model (and also models for 
polyacetyleneEHI) . While the magnetism in cuprates is, 
very well described by the SU{2) symmetric modelsE2l 
the experimental finding of stripes indicate that topolog- 
ical doping corresponds to topological defects of lower 
symmetry. Although the reason for this lowering of the 
symmetry is, most probably, dynamical in nature and is 
still not clearly understood, it gives us confidence that 
the t-Jz model is the right starting point for the descrip- 
tion of these systems. 

As long as one is concerned with the short-range 
physics of the charge and spin excitations, the isotropic 
SU{2) t-J and anisotropic t-Jz models lead tp-similar re- 
sults, as is well known from earlier studies.llj Roughly, 
the hole motion is fast and the spin relaxation is slow 
when t ^ J. Therefore, in the fast timescale the hole 
moves in the background of essentially static, staggered 
spins. We will show below that the stripe phases obtained 
numerically for the t-J^ model are victwally identical to 
those in the t-J model studied before.EII It has also been 
concluded, based on the Ginzburg-Landau functional ap- 
proach, that the antiphase shift of the antiferromagnetic 
order ppameter must originate from some short-range 
physics.Ej The rigidity of the 7r-shift of the antiferromag- 
netic phase across the domain wall in both numerical and 
experimental studies also argues for the short-range gen- 
esis of the stripes. 

In this paper we, therefore, study analytically and 
numerically the system of holes at the ADW in the 
anisotropic t-Jz model. While we restrict ourselves to 
the study of the strongly anisotropic limit of the basic 
t-J model, we believe that the results of our study are 
much more generic since the strongly correlated nature 
of the problem is preserved. 

Technically, switching off the transverse spin fluctua- 
tion has numerical and analytical advantages. The nu- 
merical advantages are twofold. First, the sizes of the sys- 
tem which can be handled by the DMRG method are con- 
siderably bigger. Second, one can think of the boundary 
conditions as a way of stabilizing ground states of differ- 
ent symmetries. In our case the choice would be between 
the state with and the state without an ADW. Then, one 
can consequently dope holes in the stripe (ADW) config- 
uration and study the evolution of the properties of the 
system as a function of doping, starting from a very di- 
lute limit. It is worth noting the boundary conditions in 
our case work as a very gentle instrument of controlling 
the symmetry of the ground state without affecting the 
wave- functions of the states themselves. The analytical 
advantage is the treatability of the problem. The analyt- 
ical part of this work largely relies on a previous study 
by two of us and Bishop, Ref. ||2l|], where we have cal- 
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culated the Green's function of the charge excitation at 
the ADW by a method which goes beyond the hmits of 
mean-field or perturbation theory. 

The general conclusion of the present study is that the 
stripe should be considered as a collective bound state of 
the holes with an ADW. In such a system the excitations 
are composite holon-spin-polarons which populate an ef- 
fective ID band. This picture is in very good agreement 
with the numerical results and provides insight into the 
problem of the stripe phase in cuprates. 

This paper is organized as follows. Section II de- 
scribes various aspects of our analytical approach in de- 
tail. Section III describes the DMRG method. Section 
le results and comparison. Section ^ lists 
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presents t 



our conclusions. 



II. ANALYTICAL APPROACH 

Our starting point is the t-Jz model which is given by: 
^ = E (45.. + H.c.) + jY.[StS^^ - -K,N,] ,(1) 



(ij) 



where t is the kinetic energy, J is the antiferromagnetic 
exchange, and Ni = n^i -I- n^. All operators are defined 
in the space without double-occupancy of the sites. 

The single-hole problem for the t-J {t-Jz) model in a 
homogeneous antiferromagnetic background is well stud- 
ied with analytical r esult s a nd numerical data being in 
very good agreementll2llH2H2a. The charge quasiparticle 
is understood as a spin pplaron, i.e. a hole dressed by 
strings of spin excitationsEa. It is also often expressed 
as that the hole movement in a homogeneous antifer- 
romagnetic background is frustrated because of the tail 
of misaligned spins following the hole, see Fig. ^ The 
idea that an ADW can be a more favorable configura- 
tion for holes relies on the fact that such a frustration 
of the hole's kinetic energy can be avoided for a move- 
ment inside the wall, such that the hole is essentially 
free in the ID structure. However, as we show below, 
the spin-polaron aspect of the physics of the charge car- 
rier remains very important in the stripe phase as well. 
We will first consider the specifics of the hole behavior 
in the inhomogeneous antiferromagnetic state (state with 
an ADW) and will address the spin-polaron aspect of the 
problem later. For the detailed description of the spin- 
polaron formalism used in this work we refer to Ref. pol . 



A. Holon in the domain wall 

Since the stripe corresponds to an ADW in the spin 
background, one has to study the nature of charge exci- 
tations at such a domain wall. Let us consider the empty 
system first. The ground state is, of course, given by 
the simple Neel configuration of spins. Fig. ^(a). How- 
ever, when the antiphase shift of the staggered magneti- 
zation is created (enforced by the boundary conditions. 
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FIG. 1: (a) A hole followed by the "string" of spin defects 
in a homogeneous AF. (b) same as (a), and "— " denote 
the sign of the staggered magnetization Mi — ( — l)'S'f . Arcs 
denote "wrong" (ferromagnetic) bonds. 



for instance) the straight bond-centered domain wall is 
the lowest energy state. Fig. |^(b). It has the energy 
^bond ~ "^/^ P^'^ ^^^^ length, which is lower, for exam- 
ple, than the corresponding energy for the site-centered 
domain wall, S/jjg — 3 J/4. We remark here that the 
SU(2) Heisenberg spins would prefer a continuous un- 
twist from one end of the crystal to the other without a 
sharp domain wall. 

Now let us consider a single hole doped to the system. 
When the kinetic energy is neglected {t = 0) the lowest 
energy state is defined from simple bond counting. Ev- 
idently, the hole is attracted to the domain wall since 
the potential energy (J-term) is lowered when the hole 
removes the "wrong" bond. Fig. ||c. 

When the kinetic energy is taken into consideration the 
following observation can be made. If one restricts the 
hole motion to one side of the ADW "ladder" of defects 
(along the y-axes, x is fixed at xq in Fig. ^ the problem 
is identical to the hole motion in the ID Ising chain.lEl 
That is, one can see that after the first step a ID spin 
defect (spinon) is created and then the motion of the 
hole does not cause any further disturbance in the spin 
background. Figs. |(a,b). The hole simply rearranges the 
"wrong" bonds while moving. Since the charge does not 
carry any "memory" about the spin of the place where 
it was created this excitation is a holon. Notably, when 
the spinon and holon are separated, they both carry a 
"kink" or "anti-kink" of the staggered magnetic order; in 
other words, they are zero-dimensional ADWs in the ID 
chain problem. Fig. ^. 

It is more visually convenient to use staggered magne- 
tization Ali — {—lyS^ instead of the on-site magnetiza- 
tion to emphasize the opposite direction of the order pa- 
rameter in the antiferromagnetic domains of spins. We, 



or 



instead of the actual 



therefore, will often use 
direction of the spin. 

Since the spinon in our case is a finite-energy excita- 
tion {E spinon — J/'^): in the t <^ J limit the hole will 
be always bound to its site of origin. One can think of 
several possibilities to avoid the spinon creation in order 
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FIG. 2: (a) Homogeneous Neel state, (b) two domains of an AF with opposite staggered magnetization separated by the 
bond-centered ADW, (c) a static hole (t — 0) is attracted to the ADW. Arcs denote "wrong" (ferromagnetic) bonds. 



to focus on the properties of the charge excitation only. 
One can (?) assume that the crystal is finite (semi-finite) 
and create a hole at one of the ends of the chain, then 
the propagation of the holon along the chain is free, {ii) 
create a pair of holes at the nearest neighbor sites and 
then consider their motion independently, in this case 
both the kink and anti-kink are carried by holons, and 
{in) start with the empty domain wall with a "wiggle", 
one half of it misplaced by one lattice spacing from the 
other half along the x-axis. This way one has an extra 
"wrong" bond in it. Fig. ||(a), which is equivalent to 
having a ID chain along the x = xq line with the sin- 
gle spinon. Then the hole creation at one of the sites 
forming the spinon is identical to the creation of a free 
holon. Figs. ||(b,c). The purpose of these manipulations 
is to show that the single hole motion along the ADW 
can be made free provided that the spin environment en- 
sures the holonic nature of the charge excitation. These 
considerations are by no means new and were discussed 
in Refs. p^ , |6^ . We suggest calling the free motion of 
the hole within the stripe a "prepared-path" motion in 
accord with the "retraceable-path" motion for the spin 
polaron. 

As it follows from the above consideration, the wave- 
function of the single holon cannot be simply written as 
a result of an action of a single annihilation operator on 
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FIG. 3: "Longitudinal" hole motion along the ADW. The 
X = Xo line can be considered as an Ising chain. Arcs denote 
"wrong" (ferromagnetic) bonds. 



some unique ground-state wave- function, since it requires 
a rearrangement of the (semi-) infinite amount of spins. 
However, if the ground-state is "prepared" , as in the case 
(m), one can still keep the formal similarity with the 
conventional single-particle creation. In other words, if 
the holon is to be created at the site i, its wave-function 
can be written as: 



II,, 



c.,.j^o,.>-(-ir^n4..n' 
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(2) 



where the Hilbert space at each site is restricted to single 
occupancy, a — —a, aj =1 (J,) if j £ A{B) sublattice, |0) 
is the vacuum state, "+" and "— " denote the domains 
with staggered magnetization "up" and "down" , respec- 
tively. This definition suffices for our consideration. 

The above mentioned "restriction" of the hole motion 
along the ADW can be formally expressed as a separation 
of the kinetic energy term in the Hamiltonian (^ in two 
pieces: 
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FIG. 4: Propagation of the hole in the Ising chain. Numbers 
indicate the amount of hoppings made by the hole away from 
its origin, "-)-" and "— " represent the sign of the staggered 
magnetization. 
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FIG. 5: (a) An empty ADW with the "wiggle". For this configuration x = xo line is an Ising chain with the single spinon. 
(b) A hole created at the place of the spinon is a holon. (c) same as (b) with "+" and " showing the sign of the staggered 
magnetization. Arcs denote "wrong" (ferromagnetic) bonds. 



where the first term ("longitudinal") includes only ID 
motion along y-axis for x — xq, and the second term 
("transverse") includes the rest. Evidently, 

nhu) = <(|i.-i) + |i.+i)) , (4) 

and thus 

nl\k) = 2t cosfc|fc) , (5) 

where 

\k)=Y,e^'^^'\U) . (6) 

i 

Consequently, the single-particle "bare" Green's function 
of the free spinless fermion (holon) propagating along the 
ADW with simple tight-binding dispersion can be written 
as 

G(°)(fc„C^) = {ky\^—r\ky) = \ ^ -n ' 

u; ^Ji'l — 2t COS ky + lO 

where the zero of energy is set at the energy of the lowest 
static {t = 0) single-hole state Eq — {H,j)t=o — Eising + 
{Ly — l)J/2 -I- 2 J, and Ly is the size of the plane in 
the direction of the stripe. The holon band minimum is 
located at ky = tt, and the index xq corresponds to the 
x-coordinate of the stripe. 

This observation that the hole can avoid the frustra- 
tion of the antiferromagnetic background in the presence 
of an ADW in comparison with the homogeneous Neel 
state where the hole motion always leads to string-like 
spin defects has been known since the discovery of the 
stripe phases. However, the sufficiency of this effect alone 
to justify the stripe formation has been exaggerated. 
Roughly speaking, this argument is gathered from the 
unphysical limit of the model, t <^ J , where the kinetic 
energy of the hole is indeed lower in the state with an 
ADW {E^^^ ^ —t) than in the homogeneous state with 



a spin polaron (i?^f„ ^ —t^/J). However, the energy 
cost of the domain wall, E^aii ^ per unit length over- 
whelms this gain in kinetic energy in this limit. More- 
over, the true ground state in the t <^ J limit is neither 
a stripe nor a homogeneous Neel state with holes, but 
rather the phase separated hole-rich and no- hole states. 
In the physical limit t ^ J the "longitudinal" kinetic 
energy of the holon at the bottom of the ADW ID band 
is —2t, while the energy of the hole (spin polaron) in the 
"bulk" (homogeneous Neel) is 

Esp = -2V3t + 0{J^/^t^^^) . (8) 

Evidently, the "unfrustrated" kinetic energy alone is in- 
sufficient to compete with the energy of the homogeneous 
state. We further illustrate this statement below. 

Figure |^ shows the energies of the spin polaron in the 
bulk and the dispersion of the pure holon state at the 
ADW for the reahstic ratio of J/t = 0.4. In the t-Jz 
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FIG. 6: The free holon band {Ek = 2tcosk, solid curve), 
energy of the spin polaron in the bulk (solid line) , and energy 
per hole in the ID holon band at half-filling, £1/2 = J/2 — 4/n 
(dashed line), are shown. 
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FIG. 7; (a) An example of a "string" generated by the "transverse" hole movement. Notably, the first defect in the "string" 
is a spinon in the x — xq Ising chain, (b) A schematic result of the hole departure from the ADW, "+" and "— " showing the 
sign of the staggered magnetization as before, (c) same as (b) with the actual spin directions shown. 



model the dispersion of spin polarons is small and we 
neglect it from this picture. One has to bear in mind 
that there is an (infinite) energy offset between these two 
Unes: the energy of the ADW Ej^^ = LyJ/2. This 
simply means that the single holon cannot compensate 
the price for the domain wall creation. Thus, the holon 
band must be "filled" to a certain level in order to reduce 
the energy. When we add more holes they will fill the 
higher /c-states in the band. Assuming a rigid-band filling 
the total energy per hole as a function of the ID filling 
fraction n|| can be calculated 

J / 1 \ 2tsin7rn|| , ^ 

Etot = - - - 1 1 . 9 

which is infinite at n|| ~f 0, Etot — for completely 
filled band n|| — 1, and has a shallow minimum at some 
intermediate value of For a chosen value of J/t = 0.4 
this minimum is around — 0.32. This lowest energy 

Emin — — 1.255i is a bit lower than the energy of the half- 
filled band E1/2 = J/2 - At/n c± -1.07t shown in Fig. | 
by the dashed line. One can see that these energies are 
more than t higher than the energy of the spin-polaron 
system in the homogeneous Neel state Egp ~ ~2.37t. 
Therefore, the energy balance of a "narrow" , strictly ID 
stripe v.s. polarons is strongly against the stripe. 



B. Transverse hole motion 

This illustration brings up the importance of the 
"transverse" part of the kinetic energy for the stripe for- 
mation. The "transverse" motion of a hole from the 
ADW, which includes all possible paths and not only 
those perpendicular to the stripe, is by no means dif- 
ferent from the "string" type of propagation in the ho- 
mogeneous AF, compare Figs. || and |3(a). That is to 
say that the charge excitation must essentially regain its 



spin-polaron properties away from the domain wall. Our 
Fig. |^(a) shows an example of a "string" generated by 
such a transverse movement. It is well known that the 
hole can propagate by erasing the taiLof "wrong" spins 
via the so-called Trugman processes.Ej There is an im- 
portant qualitative feature of our case which makes it 
different from the homogeneous problem in this aspect. 
Since the excitation inside the ADW is a holon, that is 
Q = 1, 5^ = excitation, while the spin polaron is a 
"normal" quasiparticle, Q = 1, = ±1/2, the conser- 
vation of the quantum numbers requires the departure of 
the holon from the stripe to be always accompanied by 
the emission of the spinon {Q — 0, S'^ — ±1/2). This is 
clearly the case as is shown in Figs. 0(b,c). 

In other words, the "transverse" hole motion should be 
considered as a decay process of the ID (ADW) charge 
excitation into a ID spin excitation and a "bulk" charge 
excitation, Fig. ||(a). Since both the holon dispersion 
and the holon-spin-polaron coupling are given by the 
same parameter t such virtual decays will lead to a strong 
renormalization of the holonic energy band. 

As shown in Eq. (^ the "bare" Green's function of 
the hole residing at the ADW is the Green's function 
of a free spinless fermion (holon) with tight-binding dis- 
persion. The renormalization of this Green's function 
is given by the self-energy schematically shown in Fig. 
^(b). It is well known in the single- hole problem for the 
t-J model that the self-consistent Born approximation, 
which is equivalent to the retraceable path approximation 
in our case, acccpjits for the absolute majority of such a 
renormalization.LJ This latter fact is related to the effec- 
tive analog of the Migdal theorem in this class of prob- 
lems: all first-order corrections to the hole-magnon ver- 
tex are zero because of spin conservation.LJ Corrections 
from tie higher-order processes, also known as Trugman 
paths,t3 remaining beyond the retraceable path approx- 
imation are negligible. In anticipation of the further re- 
sults we have to remark that in our problem the Trug- 
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(a) 



spinon, Q=0, S=-l/2 



holon, Q=l, S=0 




spin-polaron, Q=l, S=l/2 



(b) 




FIG. 8: (a) Decay of the holon into a spinon and a spin po- 
laron. (b) The self-energy associated with such a decay. 



man processes are simply forbidden by the conservation 
of quantum numbers for the holon decay. Therefore, the 
retraceable path approximation should work even better 
due to the fact that the holon must emit the spinon with 
the energy J. In other words, at least one spin excita- 
tion is always created and thus the Trugman paths belong 
to the same class as all other renormalization processes. 

In our case the calculation of the self-energy in Fig. 
^(b) is particularly simple since the spinon is a disper- 
sionless excitation and the double line corresponds to the 
spin-polaron Green's function known from the previous 
studiesEj. As in the case of the spin polaron, Ref. [^ , 
the renormalization is coming from the retraceable path 
movements of the hole away from ADW and back. The 
full Green's function is then given by 



1 



— 2t COS ky — Tj{u)) + iO ' 



where '^{u;) takes the form of a continued fraction. 



UJl 



UJ — UJi —UJ2 — 



(10) 



(11) 



z = 4 is the coordination number, uji is the energy of the 
I'th segment of the string, which is equal to the number of 
"wrong" bonds (J/2 each) associated with the segment, 
index xo corresponds to the a;-coordinate of the stripe. 

The energy spectrum of the elementary excitations is 
given by the poles of the Green's function Eq. ([l^) , there- 
fore one needs to calculate S(w) and seek solutions of 
E{ky) — 2t cos ky — E{E{ky)) — 0. The resulting effective 
ID band for the composite holon-spin-polaron excitation 
has been calculated in the previous work, Ref. A 
standard simplification in the calculation of S(w) is to 
assume that all oj^'s in Eq. ( pT| ) are identical so that the 
energy of the string is simply proportional to the length 
of the path and is independent of the path of a hole. 
This is a good approximation for the spin polaron be- 
cause only very few strings do not follow this rule. With 
this assumption the solution for the self-energy can be 
found in a compact analytical form given by the ratio of 



the Bessel functionsEa. If we take uji = J/2 (energy of the 
spinon) and uji^i = J [two "wrong" bonds pj er . s egment 
of the string, see Fig. 0(a)], the self-energy isE3E3 



I](a;) 



2t^ 



Lo- J/2 + V3tT{Lu - J/2) ' 



(12) 



with T(lj) = J-u./jir)/J-^/j-i{r), Mr) the Bessel 
function, and r = 2^/3t/ J. 

For our problem this "path-independent energy" ap- 
proximation also assumes that one can neglect the renor- 
malization of the "bulk" spin-polaron wave-function due 
to the vicinity of the ADW. However, there is a concern 
that this assumption may overestimate the energy of the 
excitation. Near the ADW there is less energy required 
to create a spin flip, therefore, there should be a subset 
of strings having a lower energy than the string of the 
same length in the bulk. That is, there is a question 
if the modification of the spin-polaron wave-function is 
really negligible. While we will demonstrate below the 
adequacy of the above approximation, it can be shown 
that one can consider the problem more rigorously using 
the same approach by taking into account the energy of 
each string exactly up to a certain length Ic and applying 
the path-independent assumption only for I > Ic- Tech- 
nically, it means that one can find all uji in Eq. ( pd] ) up 
to some length of the string, use the explicit continued- 
fraction form of S(a;) from Eq. ( pi] ) up to the same 
length Ic, and then use the approximate solution for the 
continued-fraction Ti{uj) from Eq. ( p2|) for I > Ic- In 
other words, the modification of the spin-polaron wave- 
function in the vicinity of the ADW can be consistently 
taken into account within the same approach. 

The results of such calculations for a representative 
value J/t = 0.4 with Ic = 4: are shown in Fig. |[ The 
energy of the lowest pole of the Green's function versus ky 
is shown by the solid line. As in Fig. ^ there is an infinite 
off-set of the renormalized holon-spin-polaron band from 
the spin-polaron energy by the magnetic energy of the 
domain wall, Et^,aii = {Ly — l)J/2. Since the effective 
band is significantly narrowed in comparison with the 
free band the energy range shown in Fig. ^ is smaller 
than that in Fig. One can see that the energy of 
the ID ADW excitation with the energy of the magnetic 
background subtracted is now lower than the energy of 
the spin polaron in the bulk at all ky which means that 
the "dressed" ID band is better for optimizing the kinetic 
energy than the homogeneous spin-polaron state. 

Another informative quantity, the residue of the 
Green's function Z{ky), is shown in Fig. || (inset). It 
gives a measure of the amount of "bare" holon in the 
wave function of the elementary excitations. We avoid 
the use of the "quasiparticle residue" since the holon is 
not a quasiparticle in a strict sense, since its overlap with 
the physical electron is zero. However, the "residue" pre- 
serves its original meaning, which refers to the renormal- 
ization of some "bare" excitation. One can see that a sig- 
nificant part of the initial holon at fc^ = tt resides inside 
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FIG. 9: Energy of arr ADW elementary excitation E{ky) vs 
ky (solid curve). Solid straight line is the energy of the spin 
polaron in the bulk E"^ . These energies are relative to the 
energy of the static hole in a corresponding magnetic back- 
ground, that is, the absolute energy of the ADW excitation 
is Email ~ {Ly — 1) J/2 higher than that of the spin polaron. 
Dashed line is the total energy per hole in the half-filled ID 
stripe band E1/2, energy of the pure Ising state being sub- 
tracted. Since the reference energy for i7i/2 and E"^ is the 
same, they can be compared directly. Inset: residue of the 
Green's function versus ky (solid curves). 



the wall (about one-half) and that almost all its weight 
is transferred to strings at jfcj, — 7r| > 7r/2. Because the 
band is very flat at the same \ky — 7r| > 7r/2, the veloc- 
ity of the elementary excitations is much slower than the 
bare Fermi velocity Wp-j^ 2t sin kp/h, in agreement with 
experimental findingsE3 It is easy to show that the veloc- 
ity at the Fermi level for our ID band is vp = Vp Z{kp). 
Generally, both the flatness of the top of the band and 
vanishing quasiparticle residue argue that the excitation 
at those k are only weakly attached to the stripe. The 
closeness of the "bulk" energy band to the top of the ID 
stripe band suggests that the ID excitations at the top 
of its band are not too different from the "bulk" spin 
polarons in the homogeneous AF. 



C. Stripe energy 

For the chosen value of J/t = 0.4 the energy of the 
holon at the bottom of the band is lower than the en- 
ergy of the spin polaron by about 1.5J, that is the en- 
ergy of three "wrong" bonds. Then the rough estima- 
tion gives that the magnetic energy of the domain wall 
will be compensated when in average every fourth site 
along the initial ID chain is taken by the hole. That is, 
the kinetic energy will make the stripe to be the ground 



state at ID linear hole concentration 



n|| > n 



1/4 (for 

J/t = 0.4). A more accurate consideration should have 
the band-filling effects taken into account. Within the 
rigid-band approximation the total energy of the system 
per hole with the energy of the pure Ising state subtracted 




FIG. 10: Total energy of the system with ADW per hole 
versus ID hole density ny (curve). Horizontal line is the total 
energy of free spin polarons in the homogeneous AF per hole. 
In both cases the energy of the pure Ising state is subtracted, 
J/t = 0.4. 



IS 



Etot/Nh 



- 1 



27m 



TV + kf 



E{ky)dky . (13) 



where the first term is the domain wall energy and the 
second term is the kinetic energy of the free "quasipar- 
ticles" filling the effective ID band up to kp = 7rn||. 
The total energy per hole of this ID band indeed crosses 
the spin-polaron energy at around n|| ~ 0.3 as shown 
in Fig. |lO|. The total energy per hole at half fill- 
ing, i?i/2 — —2.57t, is shown in Fig. ^ by the dashed 
line. Above tm. ~ 0.5 the energy (|l3| ) versus n|| is al- 
most constantEil with the energy difference between the 
stripe and no-stripe spin-polaron states of about — 0.6J. 
This behavior of Etot v.s. the absolute value of 

the energy difference, as well as the value of njj are, of 
course, functions of the J/t ratio. However, these quan- 
tities only weakly depend on J/t in the realistic range 
0.1 < J/t < 0.5 with njj shifting towards zero for smaller 
J/t (see also discussion of Etot v.s. ny in Sec. ^ ). 

An interesting question is what kind of energy scale 
defines the difference between a homogeneous state with 
spin polarons and a stripe state at some fixed n|| (say 
ny = 1/2), A£: = E'P - E^/^^. By letting J/t ^ we 
observe that AE ^ as well, and thus it cannot be pro- 
portional to t. This is in agreement with the expectation 
that at J = all magnetic configurations should become 
degenerate. It is nevertheless hard to tell analytically or 
numerically if the leading term in AE scales with J^^^ 
(however, see the discussion at the end of the section). 



D. Densities, electron distribution function 

The knowledge of the single-particle Green's function 
not only allows us to calculate the energetic stability of 



9 



the stripe as a function of n\\ and but also en- 

ables us to find spatial profiles of the charge and spin 
densities, N{x) and \ {S^){x)\, as well as tte electron dis- 
tribution function for the stripe stateE3. To evaluate 
those quantities one needs to know the wave function of 
the ID band excitation written in terms of the opera- 
tors of spins and holes. If such a wave function is known 
for each ky in the effective band then the distribution of 
the hole density and the "damage" of the spin density 
from the hole and "strings" can be extracted. Within 
the rigid-band picture these wave functions are orthogo- 
nal at different ky and thus the density profiles are given 
by the superposition of the contributions from each ky: 



N{x) = ri'''{x,ky) , 

ky<kp- 

\{S')i^)\^ E l(^'>(^'^.)l ' (14) 



ky<kF 



where 



n {x,ky) — (Jky\CjC^^\ky), 

\{s^){x,ky)\^{ky\S^^\ky) , 



(15) 



where j = {x, y) and the choice of the y-coordinate is 
arbitrary. 

Such a wave function for our composite excitation can 
be written in a variety of ways. The simplest one is to 
write it as a linear combination of the components with 




FIG. 11: The hole density modulation within the single-hole 
fcj,-eigenstate for different momenta ky, k = — 7r|. One can 
see that at the momenta away from the bottom of the band 
the hole is mostly spread around the stripe. 



strings of different length! 



^000: 



(=0 



CkyAhky) 



(16) 



where |0, ky) is a pure holonic state and under the index 
I we also understand a summation over different paths 
of the same length I. Then, using a retraceable-path 
approximation, it is easy to see that for each specific 
path 



(17) 



where multiplier \fZk'y comes from the normalization 
condition {ky\ky) — 1. An equivalent expression for the 
spin-polaron case has been obtained in Refs. |^,^ us- 
ing the spinless-fermion-Schwinger-boson representation 
for the original constrained fermion operators. 

Since the hole density and the spin density are given 
by the averages of the local (on-site) operators their cal- 
culation using Eq. (|l^) is quite straightforward because 
they are diagonal in the string basis and do not provide 
"transitions" between different components of the wave 
function ( [l^ ) . Calculation of nt is technically more cum- 
bersome since it is given by the non-local averages. One 
can always rewrite as 



where we drop the spin index for clarity. One can see 
that the matrix elements between different "strings" are 
essential. For the details of calculations of rik in the t-J 
model we refer to the earlier works, Refs. [^^8|6^. 

As an example of the calculation of the discussed quan- 
tities we show our results for n''(a;, ky), Eq. (^5|), in Fig. 
|TT| for J/t = 0.4. This quantity shows the hole density 
modulation in the direction perpendicular to the stripe 
within the fcj^-eigenstate in the ID band for different mo- 
menta ky. One can see that at the momenta away from 
the bottom of the band the hole is mostly spread around 
the stripe. Further results on these and other quanti- 
ties are given in the context of the comparison with the 
numerical data, Sec. IV. 



E. Other questions 

There are few other issues we would like to address 
here. First, if one is to start from the strictly ID stripe 
and is to gradually "relax" the initial restriction for the 
hole to be kept in such a ID chain, how soon does one 
reach the "true" eigenvalue of the t-Jz Hamiltonian? 
That is, let us allow the "bare" holon to have an ad- 
mixture of strings of length 1 = 1,2,... and see when the 
results converge and become independent of I. Figure 
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FIG. 12: Energy of an ADW elementary excitation E{ky) 
v.s. ky for "strings" with I = 0, 1,2, ...cxd included in the 
wave function. 



shows the energy bands for such "restricted-string- 
length" problems. I = corresponds to a free holon, Fig. 
^, / = 1 and / = 2 are for the cases where the "strings" of 
length 1 and 2 are included in the approximation. I = oo 
is our answer from Eqs. (|ic|)-([l2|), Fig. ^. 

This comparison demonstrates one more time that the 
"bare" , or almost "bare" ID excitation with some ad- 
mixture of the "transverse fluctuations" is certainly not 
a good approximation for the description of the stripe 
charge excitation. It is especially true for the finite con- 
centrations where large momenta around ~ 7r/2 are 
most important. 

As an alternative starting approximation we argue, 
once again, for the picture of a holon strongly coupled 
to the "bulk" spin-polaron excitation. Fig. || The dis- 
persion obtained without renormalization of the spin- 
polaron propagator in Fig. ||(b) is shown by the dashed 
line in Fig. |l^, which is already very close to the result 
where such a renowpalization is taken into account (solid 
line, from Fig. p|)llj. 




(ky-7l)/7I 

FIG. 13: ID energy band for the composite excitation. Renor- 
malization of the spin-polaron propagator at the vicinity of 
the ADW is included (solid line), not included (dashed line). 



The second issue is the following. Since there are two 
types of kinks in the system (kink and anti-kink) which 
can be associated with the pure holonic excitation. Fig. 
^, one can have two different "species" of excitations 
having opposite "geometrical" quantum numbersE3. One 
may then think of them as essentially independent par- 
ticles populating (quasi-) independent bands. 

A little thinking brings such a logic into a paradox: 
the case when the chain is completely filled (all sites are 
occupied by holes) in this language will correspond to 
two half-filled (conducting?) bands of different "flavor" . 
One resolves the paradox by noting that these "extra" 
quantum numbers should not lead to the increase of the 
total number of states in the system. In our case there 
cannot be more than one hole at any site while the as- 
sumption of "independence" of opposite species implic- 
itly doubles the Hilbert space. Therefore, if one wants to 
keep two "flavors" in the problem, an infinite on-site re- 
pulsion between them should be introduced. That latter 
problem in one dimension maps exactly onto the problem 
of spinless, single-flavor particleald and thus the single- 
band approach used throughout this work is justifled. 
Note that in the presence of the 2D degrees of freedom 
where the stripe as a whole can make wiggles the problem 
of these "gepmetrical" quantum numbers becomes more 
complicatedlS. 

The last question is what energy scale deflnes the 
"free" kinetic energy difference between the stripe and 
no-stripe states. Naively, such AE should be governed 
by the kinetic t-tevm since the hole motion can be made 
free in a ID structure. However, the energy of the or- 
der ~ JL is paid to "prepare" such a structure in a 
2D AF. Let us consider the t ^ J limit for the ID 
hole motion and assume that the length of the ADW 
is a free parameter which would minimize the total en- 
ergy. In the continuum limit {Ekin) = —2t-\- At/L'^, 
(Ej) ^ JL, and the minimum of the energy is achieved 
at Lopumai ^ {J/ty^^. The corresponding energy is 
Emin — — 2i -I- a{J^ty/^. One recalls an almost iden- 
tical consideration of the "retraceable-path" motion of 
the hole by the "strings" in the spin-polaron problem 
which also gives (Lstring) ~ {J/tV^"^ and Esp ~ 2^/3t + 
l3{j'^t)^^^, see Ref. Therefore, there is no new en- 

ergy scale, different from the homogeneous Neel problem, 
introduced by the domain wall and the "prepared-path" 
motion in ID ADW is, in fact, not too different from the 
"retraceable-path" motion of the hole within the spin po- 
laron. One may conclude that the same scale ~ (>/^i)^A 
should govern the energetic balance favoring the stripe.O 
This is yet another argument that the stripes are the out- 
come of the same tendencies which are seen already for 
the single-hole problem. 



III. NUMERICAL APPROACH 

The numerical study of the t-J^ model utilizes the 
DMRG method for the clusters up to 11 x 8 with var- 
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FIG. 14: 11 X 6 cluster with one hole, cylindrical BC's, and 
J/t — 0.35. The size of the circles is proportional to the 
on-site hole density. 





FIG. 16: 11 X 4 cluster, one hole, open BC's in the narrow 
direction, no staggered field, cylindrical BC's, J/t — 0.2. The 
ground state here contains a stripe. 




FIG. 15: 11 X 6 cluster with one hole, field applied at open 
ends (crosses) to stabilize ADW, J/t = 0.35. 



FIG. 17: 11x8 cluster with four holes, cylindrical BC's, and 
J/t — 0.35. The ground state contains a stripe. 



ious boundary conditions (BC's). The BC's at the left 
and right sides are always open; if the BC's at the top and 
bottom are periodic we refer to the BC's as cyhndrical. A 
staggered magnetic field, typically of size O.lt, can be ap- 
plied at the open ends in order to enforce the ADW inside 
the system, or enforce it to stay homogeneous. Typically 
we keep 1000-2000 states per block, and perform up to 
a dozen finite system sweeps. Typically the truncation 
error was 5 x 10~^. 

Figures |l^0 show several exainples of the states 
within the different clusters. Figure |lj shows the 11x6 
cluster with one hole, cylindrical BC's, and with J/t — 
0.35. The ground state is homogeneous in this case. 
Figure |l^ shows the same 11x6 cluster with one hole, 
but with a staggered field applied at the open ends, 
J/t = 0.35. The ground state is an ADW with the hole 
bound to it. Figure [l^ shows the narrower 11 x 4 cluster 
with one hole, no staggered field, cylindrical BC's, and 
at smaller J/t = 0.2. The system spontaneously forms 
a stripe in the ground state for this system. One may 



conclude from here that in the 11 x Ly system with the 
number of holes Nh ~ Ly/A one will have a stripe in the 
ground state for J/t < 0.2. Figure |l^ shows the largest 
11x8 cluster with four holes, no staggered field, cylin- 
drical BC's, and J/t = 0.35. In this case the ADW is 
also spontaneously formed with the holes occupying the 
border between the antiferromagnetic domains. 



IV. RESULTS 

A. Single-hole problem 

The first problem we would like to clarify is the single- 
hole excitation problem in the stripe configuration. We 
therefore consider different clusters with an ADW in- 
duced by staggered fields for various values of J/t. 
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FIG. 18: A schematic picture of a 7 x 5 cluster with periodic BC's and ADW. (a) One hole, straight ADW is enforced by 
periodic BC's, (b) holon and spinon in the straight ADW, (c) spinon and holon recombined after the circular motion of the 
hole around the cylinder, the ADW is translated in a;-direction as a result. The shaded row at the bottom is equivalent to the 
top row and is shown to emphasize periodic BC's in y-direction. 



1. Boundary conditions 

We would like to remark here that simple cylindrical 
or periodic BC's do not quite work for the purpose of the 
study of the single excitation for the following reason. As 
we have discussed in Sec. II the "core" of the composite 
excitation is a holon, that is, an excitation carrying the 
kink. Fig. ||, while periodic BC's require the domain wall 
to be closed on itself around the cyHnder, Fig. p^(a). 
This means that (i) there will be two excitations in the 
wall, holon and spinon, instead of one. Fig. |l^(b), (ii) 
the free holon movement is frustrated by the spinon, [Hi) 
periodic BC's induce an artificial (absent in — > oo 
limit) stripe meandering. The meandering is induced by 
the process of holon recombination with the spinon with 



the subsequent change of the side of the domain wall. 
Then, when the holon completes the full circle around 
the cylinder the whole domain wall is translated in the 
x-direction. Fig. |l8|(c). 

The way to avoid such a frustration of the "free" 
holoniC|— motion in the finite system is intuitively 
evidentO. Since the holon carries a topological ("geo- 
metrical") charge the topological BC's are to be used. In 
our case it is Mobius BC's, Fig. ^ One can see that 
with these BC's the holon motion is unfrustrated. The 
number of sites should be odd to avoid the frustration 
in the spin system. For the case of the Mobius BC we, 
therefore, study clusters 11x7,11x5, and 11x3, and 
in the case of cylindrical BC's, the clusters are 11 x 8, 
11 X 6, and 11 x 4. 




FIG. 19: A schematic view of 7 x 7 cluster with Mobius BC's. 
An ADW in the center of the cluster and a "wiggle" are in- 
duced by Mobius BC's. 



2. Energy 

Figure ^ shows the J/t dependence of the ground- 
state energy of a single hole in the 11 x 7 cluster with a 
Mobius BC, energy of the undoped system being sub- 
tracted. The theoretical curve shows the J/t depen- 
dence of the ground-state energy of the single holon-spin- 
polaron, that is, the energy of the bottom of the effective 
ID band in Fig. |, obtained from Eqs. (0)-(©. Note 
that the reference energy in Sec. II differs from that used 
here by the constant 2 J, the quantity associated with the 
energy of four bonds "broken" by the hole. We would like 
to note that there is no energy adjustment between nu- 
merical and analytical data used in Fig. 20 and it is 
not the result of the best fit. The maximal discrepancy 
of numerical and analytical results in Fig. ^ is about 
~ 0.1%. The closeness of agreement is even better than 
in the case of similar calculations for the ^n polaron 
compared to the Exact Diagonalization dataE3. The rea- 
son for that has been discussed in Sec. II. The Trugman 
paths, which are beyond the retraceable path approxi- 
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FIG. 20: J/t dependence of the ground-state energy of: the 
single hole in 11 x 7 cluster with Mobius BC (DMRG, circles) 
and the single holon-spin-polaron in the infinite stripe (theory, 
solid line). 



mation, are the source of the discrepancies in the case of 
the spin polaron. In our case the Trugman paths do not 
exist in their original form. 

One can also study the ground-state energy depen- 
dence on the size of the system. In Fig. |2l| we show 
Egs versus 1/Ly for Ly = 3,5,7 (Mobius BC's) and for 
Ly = 4,6,8 (periodic BC's) together with the theoreti- 
cal results for the infinite system, J/t = 0.35. For the 
case of Mobius BC's the theoretical point at 1/Ly — 
is Egs — Ehoiom for the case of periodic BC's the theo- 
retical point represents the sum of the holon and spinon 
energy, Eqs = Ehoion + J /'2- That is, we assume that 
they do not form a bound state in the thermodynamic 
limit. One can see essentially negligible finite-size effects 
on the energy of the "free" holon-spin-polaron (Mobius 
BC's, also see inset). This feature can be anticipated 
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FIG. 21: The ground-state energy v.s. inverse linear size of 
the cluster. Mobius BC's (circles) correspond to the case of a 
single holon, periodic BC's (squares) induce the state with the 
holon and spinon. Theoretical results for the energy of the free 
holon and holon 4- spinon are put on the j/-axis (diamonds). 
Lines are guides to the eye. Inset shows Mobius BC's and 
theoretical data in a strongly magnified scale. J/t = 0.35. 
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FIG. 22: Hole density distribution across the stripe at J/t = 
0.35. Data for Mobius BC's, 11 x 7 cluster (circles), peri- 
odic BC's 11x8 cluster, and theoretical results (crosses) are 
shown. Lines are guides to the eye. Inset shows Mobius BC's 
and theoretical data for a:: = 6 in a strongly magnified scale. 
Enlarged area is shown by the box. 



since the excitation is a band-like state and the ground- 
state wave-vector belongs to the reciprocal space of all 
three clusters. The linear extrapolation to 1/Lj, = in 
Fig. ^ was made from Ly = 5,7 and Ly ~ 6,8 for 
Mobius BC and periodic BC results, respectively. 



3. Density 

Figure |2| shows the results for the hole density distri- 
bution across the stripe for the ground state of the single 
hole at J/t = 0.35. DMRG data from the 11 x 7 cluster 
with Mobius BC's are shown by the circles. Theoretical 
results obtained from Eqs. (p^)-(p^) are shown by crosses. 
The total disagreement of the numerical Mobius BC data 
and analytical result is Y.x l^«"m - Nth\/Nnur,i ^ 0.3%. 
The same Fig. ^ shows the density profile of the hole 
distribution along the x-axis for the case of periodic BC's. 
We note that the hole density modulation in the case of 
periodic BC's is essentially induced by the repulsion of 
the open boundaries in the x-direction. If the system 
would be significantly wider in the a;-direction the circu- 
lar meandering effect, which induces the transverse mo- 
tion of the domain wall, would spread the density homo- 
geneously. The same would be true in the case of periodic 
BC's in the a;-direction. This just demonstrates the fol- 
lowing paradox: a finite system with periodic BC's in all 
directions would show a homogeneous ground state, as in 
Ref. [ p^ , whose wave function may be, in fact, given by 
the superposition of slowly meandering stripes. In that 
case the stripe should be detected not from the density 
profile but from the instantaneous spin-spin correlation 
function, which should exhibit a strong anti-phase com- 
ponent. 

The density profile study of the single-hole problem is 
completed by the |(S'^)(a;)| data shown in Fig. Note 
that the spin-density profile is not straightforwardly re- 



14 




10 11 



FIG. 23: The spin-density 1(3") {x)\ v.s. x at J/t = 0.35. 
Data for Mobius BC's, 11 x 7 cluster (circles), and theoretical 
results (crosses) are shown. Lines are guides to the eye. Inset 
shows Mobius BC's and theoretical data for x — 4 and 5 in a 
magnified scale. Enlarged area is shown by the box. 



0.52 
0.51 
0.50 

j 

['0.49 
0.48 
0.47 
0.46 



-J/t=0.35, 1 hole, MBC 


1 









(OX) 


^ (k,k) 




I.I.I. 


1 



0.5 



1.0 



1.5 

k 



2.0 



2.5 



3.0 



FIG. 24: Jik for one hole in Mobius ECU x 7 cluster for (1, 0), 
(0, 1), and (1, 1) directions, J/t = 0.35. Solid and dashed lines 
are DMRG and theory results, respectively. 



4- Electron distribution function 



lated to the hole density since there is also a large contri- 
bution from the hole-induced strings of misaligned spins 
affecting the average on-site spin values. For a single hole 
in the stripe of the length Ly the suppression of the spin- 
density should be proportional to 1/Ly [Ly = 7 in the 
case of the 11x7 cluster). This is valid for all x except 
the "center" of the stripe where the motion is holonic. 
Since the holon is a topological excitation it borders two 
ID domains with opposite staggered spin value. Because 
of that any site along x = xq [xq = 6 in Fig. ^3|) al- 
ways has average spin zero and this is not a, 1/Ly effect. 
The deviation of the theoretical results from DMRG for 
|(S'^)(a;)| is more observable than for N{x). This is be- 
cause of the subset of the paths which correspond to the 
processes of hole departing and then rejoining the stripe 
in the other leg of the ladder. Such paths can produce de- 
viations from the 1/Ly character of the spin-density mod- 
ulation for X xq, similar to that at x — xq, which are 
hard to account for analytically. In any case, the largest 
discrepancy between the theory and numerical data for 
ISS'^ix ^xo± l)/S''{x)\ does not exceed 10%. 



TABLE I: E.(4e»+d> 
J/t = 0.35. 



for several coordination vectors d. 



d 


DMRG 


theory 


(1,0) 


0.1041 


0.102 


(0,1) 


0.1216 


0.115 


(1,1) 


-0.1266 


-0.093 


(2,0) 


-0.0030 


~ 


(0,2) 


-0.0041 


~ 



The electron distribution function is compared as well. 
For this quantity the analytical calculations are quite 
cumbersome and were restricted by |d| = 2 in Eq. (|l|). 
One can see from Eq. ( p^ that, apart from the constant, 
nk is given by the set of cos(dk) with the coefficients 
given by the ^i{clci+d) averages. We fist such coeffi- 
cients for several d's in Table I, where one can see a very 
good agreement of DMRG with the theory. 

Figure ^ shows a plot of nk in the (1, 0), (0, 1), and 
(1,1) directions. A more informative 2D intensity plot 
in Fig. |2^ shows the DMRG results for the electron dis- 
tribution function averaged over the x and y directions, 
n\i = [n(kx:,ky) + n{ky, kx)]/2. Two features are worth 
discussing. 

First, there is almost no anisotropy between the k^ 
and ky directions in for the highly anisotropic stripe 




FIG. 25: 2D intensity plot of Tik ~ [n{kx,ky) + n{ky,kx)]/2 
from DMRG for a single hole in Mobius BC 11 x 7 cluster, 
J/t = 0.35. 
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FIG. 26: A schematic picture of 7 x 5 cluster with two holes in an ADW. (a)-(d) hole motion in the periodic direction combined 
with the transverse motion of both holes results in the translation of ADW in a;-direction. The shaded row at the bottom is 
equivalent to the top row and is shown to emphasize periodic BC's in y-direction. 



configuration, Fig. |2J. That can be, broadly speaking, 
interpreted as yet another demonstration of the equal im- 
portance-pf the "transverse" and "longitudinal" kinetic 
energiesO. The broad features in in the i- J-like mod- 
els are understood as coming from the "fast" (incoherent) 
motion of the hole inside the quasiparticle, while the co- 
herent part should show itself as a (J-peaL-ai-k = ko 
proportional to the quasiparticle residuet3l£3'LJ. In the 
stripe case, excitations are not quasiparticles in a stan- 
dard sense and, therefore, no sharp features are to be 
expected. In fact, it is easy to show that rtk for a holon 
in a periodic ID Ising chain is equal to: 



J^k = {ko\cl.ci^\ko) 



- — Y cosfccosfco , (19) 



where fco is the momentum of the holon. This is the 
simple consequence of the fact that the holon is a zero- 
dimensional domain wall and that (c^Ci+d) averages for 
distances larger than d = 1 are identically zero. 

Second, the general behavior of rik in our problem is 
similar to that for the spin polaron: it has a maximum 
at k = (0, 0) and a minimum at k = (tt, tt) which is a 
simple consequence |-Ojf the kinetic energy minimization 
in the ground stateS. However, its shape is cross-like 
rather than diamond-like (for the shape of spin-polaron 
rik see Refs. |^,|9j). For spin polarons nk is mainly 
given by the powers of 7k = (cos + cos fcj^)/2, which is 
zero along the lines (tt, 0)-(0, tt). This gives a diamond- 
like shape of rik- In the ADW configuration {c\ci+d) for 
d = (2, 0) and (0, 2) are strongly suppressed due to the tt- 
shift of the antiferromagnetic order parameter across the 
wall and across the holon. As a result, there are signifi- 
cant and almost equal cos k^^y and cosfc^; cos ky harmon- 
ics of opposite sign in rik (see Table I), but the harmonics 
cos 2kx^y are much suppressed. This is different from the 
spin-polaron case where they are all of the same order. 
Our results for rik in Fig. g5| are obtained for the single- 
hole problem in a stripe configuration, but since rik is 
"saturated" on short distances the same characteristic 
features should remain the same for higher doping. In 



fact, our "Maltese-cross" -like shape of rik is remarkably 
close to that observed in angle-resolved photoemission 
spectroscopy, Ref. |75| . 

Altogether, we have a very close agreement of the the- 
ory and DMRG data on the energy, density, and electron 
distribution function for the single excitation in the stripe 
configuration. This proves that our description of such 
an excitation is correct. 



B. Many- hole problem 

The second problem we address in this work is the 
many-hole system. 



Ill the fi- 



As we described in Sec. 
nite cluster can be doped with different amount of holes 
and one can use staggered-field BC's at the open bound- 
aries to enforce the state with and without an ADW. 
At some doping concentration the stripe state with an 
ADW becomes a ground-state and no fields are neces- 
sary to stabilize it. We also note that for two and more 
holes the Mobius BC's are not beneficial anymore be- 
cause the holes play the role of boundaries for each other 
and the meandering of the domain wall can come as a 
result of some "collective" motion. Our Fig. ^ shows 
an example of such a process for the case of two holes. 
Such effects were discussed earlier in Refs. 



5 63 . How- 



ever, we will show that the role of such processes in the 
stripe energy is negligible. This can be anticipated since 
the holes are spread significantly within the individual 
excitation and such collective processes should be statis- 
tically rare. Moreover, the "bending" of the stripe affects 
the free longitudinal motion of the holes and thus is unfa- 
vorable. This is in accord with the earlier work, Ref. 
where such a rigidity was referred to as a "garden hose" 
effect. Note that in the case of a strictly ID stripe close 
to complete filling (n|| ~ 1) the longitudinal kinetic en- 
ergy is suppressed, and the stripe meandering can become 
more importantE3. 
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FIG. 27: Total energy of the system with an ADW per hole 
versus n||. Circles are the DMRG results from the 11 x 8 
cluster. Solid curve and crosses are the theoretical results as 
described in the text. Horizontal solid and dashed lines are 
the energies of free spin polarons and bound states of spin 
polarons in the homogeneous AF, respectively. J/t — 0.35. 



1. Total energy 

In Fig. ^ we show the DMRG and theory results for 
the total energy of the system per hole versus linear hole 
concentration n|| = Nh/Ly (energy of the empty system 
being subtracted). DMRG data are obtained in the 11x8 
cluster with Nh — 1,2,3,4,5 and 6 holes populating an 
ADW. At Nfi > 3 there was no need to stabilize the stripe 
by the fields at the boundaries since it was the ground 
state of the system. The theoretical curve is equivalent 
to the results shown in Fig. |l^, which are obtained from 
the rigid-band filling of the effective ID band using Eq. 
(|l3|). Note again that the reference energy is different 
from Sec. II by 2J. Crosses are also the result for the 
same rigid-band filling but for the system with the dis- 
crete k-space. We mimic the periodic finite-size system 
by imposing that only Ly = 8 k-points are available in 
the effective ID band. Straight solid and dashed lines 
are the energies per hole of the systems of independent 
spin polarons and bound states of spin polarons in the 
homogeneous AF, respectively. 

Obviously, the rigid band filling neglects the effects 
of interaction between the carriers except the Fermi re- 
pulsion. Such interactions are quite complicated and 
would include attractive as well as repulsive terms as 
well as some collective effects resulting in stripe mean- 
dering. However, given the good agreement of analytical 
and numerical results one can conclude that such effects 
are secondary for the stripe formation. Therefore, the ki- 
netic energy of the individual holes, both along the stripe 
(holonic motion) and perpendicular to it (spin-polaron 
part), is the main reason that brings the stripe to the 
ground state. Note that since the holons are spinless, 
the Fermi exclusion is much more effective at inducing 
proper correlations between holons, acting as a hard core 
repulsion, than in the usual spin-1/2 case, where up and 
down spin particles can be on the same site. Thus, the 



effectiveness of this simple band filling approach is not so 
surprising. 

One can question the physical picture of a "straight" , 
weakly meandering stripe described by an effective ID 
band from the point of view of applicability of the "free- 
holon" approximation, Ref. . The controversy is that 
at the physically relevant concentration n|| ~ 0.5 the free- 
holon approximation (ny <C 1) and the free-electron ap- 
proximation (1 — 7i|| <C 1) are not applicable and should 
be equally bad. A very good agreement of our theory 
with the numerical data up to n\\ ~ 0.5 and beyond can 
be seen as quite puzzling since the theory is based on 
n|[ ^ 1 approach. However, such a controversy comes 
from the mean-field picture of the strictly ID stripe. 
One should rather consider a stripe to be a combina- 
tion of strongly dressed, well spread holes forming a col- 
lective bound state with the ADW. In fact, the actual 
amount of holons within the effective ID band is given 
by rihoions = J2ky<kF {Zky is a residue of the Green's 
function), which, for the physical range of parameters, 
does not exceed Uhoions — 0.2 <C 1 even for a completely 
filled stripe, rty — 1. It demonstrates that the rty <C 1 
approach should work well in all ranges of doping. 



2. Chemical potential 

Figure ^ provided a comparison of the total energy of 
the system within the different topological sectors, which 
defines the ground state. However, a more discrete ener- 
getic analysis is necessary to study the delicate balance 
of the stripe formation. Here we introduce the chemi- 
cal potential as the difference between the energy of the 
system with Nh and Nh — 1 number of holes: 



l^{Nh) = E, 



Nh-l 



(20) 



In the situation when the kinetic energy is frustrated, 
it gives a measure of how effectively the energy of the 
system is lowered by an extra hole for the states with 
different topology. 

Figure 28 shows /i as a function of Nh for the DMRG in 
the 11x8 cluster with cylindrical EC's for the states with 
and without the stripe, together with theoretical results 
calculated from the discrete-A: rigid-band filling of the 
effective ID stripe band. In obtaining theoretical data 
points we needed to account for the frustrating character 
of the cylindrical EC's (discussed above for one hole) for 
the cases when the number of holes is odd: 



^tot 



J2 Ek^n{kn) + (3J/2)iV,, + {J/2)d_ 



Nh .odd 



(21) 



where fc„ is one of the available Ly = 8 k-points, counted 
from the bottom of the band, Ekm{kn) is the energy of 
the ID excitation. Fig. ^, 3 J/2 is the energy of the static 



hole at the ADW, J/2 = E, 



is the frustration energy 



caused by EC's. As a result, the theoretical expression 
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FIG. 28: Chemical potential vs A^^. Discrete-fe rigid-band 
results, Ly = 8, (diamonds), DMRG data in 11 x 8 cluster for 
the system with ADW (circles) and without ADW (squares) , 
J/t — 0.35. Lines are guides to the eye. 

for the chemical potential is given by: 

fi^'iNh) = Ek^nikNj + 3J/2 + (J/2)(-l)^" , (22) 

where EkmikN^) is the lowest energy available for the 
Nh hole in the ID band. The last term provides a zigzag 
behavior of fi{Nh) shown in Fig. ^ 

For Nh — 1 the DMRG and the theory points are 
equivalent to those in Fig. ^ for the Ly = 8 cluster 
with periodic EC's and for the E^f = Ehoion + Espinon, 
respectively. Since the DMRG data in Fig. 12^ seem to 
scale to Eg^, the difference ^^^'"'^{1) - yF^'"^y{l) is, 
most probably, a finite-size effect. One can see a very 
good overall agreement of the trends in both the numer- 
ical and analytical result. Since the theoretical results 
are calculated from the picture which corresponds to the 
subsequent filling of the ID band, the chemical potential 
should grow as the higher fc-states are filled. 

We mark differently the data for Nh = 5 and Nh — 6 
since in the Ly — % stripe they correspond to the 
high concentrations where the finite-size effects become 
more pronounced. Note again that the ADW configu- 
ration (stripe) is enforced by the boundary conditions 
&i Nh = 1,2 while for Nh = 4,5,6 it is the ground 
state. Conversely, the no-stripe state (homogeneous AF) 
is the ground state for Nh = 1,2 {Nh = 3 looks like a 
metastable state) and is to be enforced for Nh > 2. We 
mark the no-stripe data for Nh = 3, 4 as shaded because 
the Nh — 4 state is not formed by the gas of spin po- 
larons or polaron pairs but rather is a stripe-like circle 
with an ADW in the center of the cluster (Fig. ^) . It is 
not a bound state of polaron pairs, but rather a many- 
particle bound state with the condensate of magnons, 
the "droplet" of the 7r-shifted AF inside the circle cor- 
responds to such a condensate. This shows the strong 
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FIG. 29: 11 X 6 cluster with four holes. "Homogeneous" an- 
tiferromagnetic state is enforced by the boundary conditions. 
Instead of being homogeneous, the circular ADW is created. 



tendency to the stripe formation such that even a small, 
finite number of holes will prefer to form a closed loop 
of the ADW "nuclei" , which can then develop into the 
straight stripes as the doping grows. We would like to 
discuss here the Nh ~ I and Nh ~ 2 cases in Fig. ^ 
in more detail. While the no-stripe state is the ground 
state for these hole concentrations the effectiveness of the 
energy lowering is much higher in the stripe state. For 
Nh = 1 it is yet another form of the discussion given in 
Sec. II that at the bottom of the ID holon-spin-polaron 
band the energy is significantly lower than in the bulk. 
In other words, from Fig. ^ one can see that the indi- 
vidual charge carries benefit energetically from being at 
the stripe. Since for Nh = 1 in either the stripe or the 
no-stripe system it is the kinetic energy of the individ- 
ual charge excitation which is optimized, one can con- 
clude that the formation of the stripe is "kinetic-energy 
driven" . 

The idea that stripe formation can be viewed as a con- 
densation of a set of hole pairs has been discussed in 
Ref. [Q. Here we see that if one preconfigures the ADW, 
even a single hole is bound to it with substantial bind- 
ing energy, and pairing is not involved. The quantitative 
description of stripe formation that we have developed, 
which does not involve pairing, suggests that pairing is 
a lower energy phenomena which can be considered af- 
ter the stripe is formed. However, the condensation of 
pairs idea may also be a valid point of view. An ener- 
getic test for this point of view would be that the energy 
per hole of a stripe should be only slightly less than that 
of separate pairs. More specifically, the energy advan- 
tage of the stripe should be less than the pair binding 
energy. For the 11 x 8, J/t = 0.35 system, we find that 
the energy of the half-filled stripe, per pair of holes, is 
2El*l'PyNh = -3.91. This compares with E2 = -3.72 
for a single pair. The difference between these two ener- 
gies is approximately twice the binding energy of a pair 
of holes, A ~ —0.09, suggesting that the condensation 
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of hole pairs is not nearly as good a description as the 
current formalism, at least in the t-J^ model. 

A separate issue is whether the stripes help to promote 
pairing, are irrelevant to it, or hinder it. In the no-stripe 
state the chemical potential is lower for Nh — 2 than 
for Nh = 1. This energy is the true bound-state energy 
of the pair of spin polarons, ^sv = — /^(l) < 0, 
much studied in the past, Ref. |20|. In the stripe state 
the excitations do not form a true bound state, although 
n{2) — fi{l) < 0. The theoretical result for this difference, 
/^(2)-^(l) ^ Ekinih) - Ekinih) - J, contains the neg- 
ative energy — J provided by the removal of the spinon, 
which is induced by the PBCs in the one-hole system and 
has nothing to do with the pairing. Since the DMRG data 
look very much the same and also ^^*^^<^(2)-/i*''^°'^?'(2) 
is almost the same as for the one-hole case, one may con- 
clude that there is no binding involved here at all. How- 
ever, it may also be that the binding energy is compen- 
sating the stronger finite-size effects for the N/^ = 2 case 



leading to the same ^,^^'^^^{2) - 



theory 



(2). In any case. 



we see no significant enhancement of binding and thus 
the stripes seem to be largely irrelevant to pair-binding. 

The collective stripe fluctuations (as opposed to the 
fluctuations of individual holes) have been discussed as 
an effective alternative-way to lower the energy and sta- 
bilize the stripe phaseEj. We check the effectiveness of 
such processes by studying two holes in the 11 x 8 sys- 
tem with periodic EC's and Mobius BCs. The Mobius 
EC's suppress the stripe meandering considerably, which 
is seen in the hole density profiles, while the difference 
between the energies of the states is slightly above the 
numerical accuracy of the DMRG. 

The overall conclusion of this discussion is the follow- 
ing. The kinetic energy of the individual charge excita- 
tions at the domain wall, which includes significant com- 
ponents of both "longitudinal" and "transverse" motion, 
is the reason for the stripe formation. The pairing en- 
ergy does not seem to be significantly modified and, in 
general, is associated with the smaller energy scale. Me- 
andering of the stripe, while being important for some 
other properties, has very little effect on the energy of 
the stripe. Therefore, one can write this as a hierarchy 
of the energy scales: 



(23) 



where AEkm is the difference in the "release" of the ki- 
netic energies between the stripe and homogeneous state, 
AEkin ~ {,Pty^^, as discussed in Sec. II. Epairing J 
but, in fact, is only a fraction of J, Ref. Emeandering 
should carry a statistically small factor describing the 
probability of the collective motion of two or more holes 
together. As we mentioned, Emeandering is hardly de- 
tectable numerically. 

The emerging picture of the stripe as a collective bound 
state of strongly dressed ID band excitations with the 
ADW also implies that there are "deep" states, which 
reduce the energy of the stripe, and "shallow" states, 
which are spread around the stripe and are only weakly 
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FIG. 30; Hole density distribution across the stripe at J/t = 
0.35. DMRG data for four holes in 11 x 8 cluster, cylindrical 
BG's (filled circles). Theoretical results for (i) a half-filled 
stripe centered at a;o = 6 (empty circles) ; (ii) linear combina- 
tion of stripes at xq = 5, xq — 6, and xq = 7 (crosses) are 
shown. Lines are guides to the eye. 



coupled to it. 



3. Density 



To conclude this section we show the hole density pro- 
file for the half- filled stripe in the 11 x 8 cluster (four 
holes) compared to the theoretical results for n|| = 1/2 
from Eqs. (]l3)-(|l|), J/t = 0.35, Fig. |o[ In DMRG data 
cylindrical EC's are used and the stripe is the ground 
state for this system. Theoretical results for the ID-band 
stripe centered in the middle of the system at a:o = 6 are 
shown by the empty circles. As we discussed before, the 
meandering effect causes the migration of the stripe as 
a whole in the transverse direction. Since we know that 
the meandering effect is weak and that it effectively leads 
to the coupling of stripes centered at the different xq we, 
therefore, model this effect by assuming that the ground 
state is given by the linear superposition of ID bands 
centered at a; = 5, 6, and 7 with equal weight. Further 
distribution is assumed to be unfavorable because of the 
open EC's. Then the density profile is given by the aver- 
age: N{x) = [N,,=5{x) + N,^=e{x) + N^^^7{x)]/3. The 
results of such an averaging are also shown in Fig. ^ 
(crosses). Even better agreement can be reached assum- 
ing wider meandering and weight distribution, but such 
a task is beyond the scope of this work. 



V. CONCLUSIONS 

Summarizing, we have presented a comprehensive com- 
parison of the DMRG numerical data for clusters up 
to 11 X 8 with the analytical studies based on the self- 
consistent Green's function method for a single stripe of 
holes in an AF described by the t-J^ model. We consider 
the close agreement of the results as a strong support for 
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the validity of our analytical method and of the phys- 
ical picture which follows from it. We have provided 
a description of the charge carriers building the stripe 
as a system of ID elementary excitations, unifying the 
features of holons and antiferromagnetic spin polarons. 
Then the stripe should be seen as an effective ID band 
partly filled with these elementary excitations. This pic- 
ture is in a very good accord with the numerical data. 

As it follows from our study, the stripe can be roughly 
described by the deep "backbone" states, which mini- 
mize the energy of the anti-phase configuration in the AF 
and the shallow, almost free spin-polaron-like excitations 
around the ADW. Since the spin polarons are known to 
have a considerable pairing between themselves, such a 
framework does not require the superconducting pairing 
to come from some ID instability, but rather suggests 
that the pairing is largely unrelated to the ID stripe pat- 
tern. Such a scenario is also discussed in recent work, 
Ref . 1?^ . Another hypothetical advantage of our picture 
is a more effective screening of the long-range component 
of the Coulomb repulsion, which represents the problem 
for the system of strictly ID charges.Ea 



Altogether, the comprehensive comparison of the re- 
sults of the theory and DMRG numerical approach has 
shown a very close quantitative agreement, thus provid- 
ing a strong support to our way of understanding the 
charge excitations at the anti-phase stripe in an AF. 
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